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Abstract. We demonstrate that the Auger decay rate in an endohedral atom is very 
sensitive to the atom's location in the fullerene cage. Two additional decay channels 
appear in an endohedral system: (a) the channel due to the change in the electric 
field at the atom caused by dynamic polarization of the fullerene electron shell by 
the Coulomb field of the vacancy, (b) the channel within which the released energy 
is transferred to the fullerene electron via the Coulomb interaction. The relative 
magnitudes of the correction terms are dependent not only on the position of the 
doped atom but also on the transition energy u>. Additional enhancement of the 
decay rate appears for transitions whose energies are in the vicinity of the fullerene 
surface plasmons energies of high multipolarity. It is demonstrated that in many 
cases the additional channels can dominate over the direct Auger decay resulting in 
pronounced broadening of the atomic emission lines. The case study, carried out 
for Sc 2+ @Cg^, shows that narrow autoionizing resonances in an isolated Sc 2+ within 
the range ui = 30 . . . 45 eV are dramatically broadened if the ion is located strongly 
off-the-center. Using the developed model we carry out quantitative analysis of the 
photoionization spectrum for the endohedral complex Sc3N@Cso and demonstrate 
that the additional channels are partly responsible for the strong modification of the 
photoionization spectrum profile detected experimentally in [1]. 
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1. Introduction 

In this paper, we demonstrate that the rate of atomic non-radiative decay in an 
endohedral system A@Cn strongly depends on the location of the atom inside the 
fullerene cage. In particular, the width of the Auger decay can be varied by orders 
of magnitudes by increasing the atomic displacement a from the cage center. As a 
increases the most pronounced enhancement of the width occurs for the transitions 
whose energy lies in the vicinity of the fullerene surface plasmon of high multipolarity. 
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To analyze the influence of the fullerene electron shell on the atomic decay we 
consider two physical mechanisms. The first one accounts for the change in the electric 
field at the atom due to the dynamic polarization of the shell by the Coulomb field of 
the transferring atomic electron [2]. This scheme is similar to the one which leads to 
the modification of radiative and Auger decays in an isolated atom due to the intrashell 
many-electron correlations [3] . Within the second mechanism the energy, released in the 
atomic transition, is transferred via the Coulomb interaction to the fullerene electron 
which becomes ionized. In [4] this mechanism was called interatomic Coulombic decay, 
thus pointing out its close relationship with the interatomic decay process in various 
molecular clusters [5-7]. 

The influence of the fullerene electron shell on the processes of radiative and non- 
radiative decays of a vacancy in the atom was analyzed for spherically symmetric 
endohedral systems A@Ceo, in which the atom A is placed at the center of the spherical 
cage of the C 60 fullerene [2,4]. It was noted that the interatomic Coulombic decay 
leads to a dramatic increase of the atomic decay rate, especially in the case when a 
non-radiative transition in the isolated atom is energetically forbidden. The estimates 
presented in the cited papers show that the interatomic decay mechanism in A@Ceo can 
increase the radiative decay rate by a factor 10 5 — 10 6 . On the other hand, the impact 
of the dynamic polarization of the fullerene on the atomic Auger width was found to be 
negligibly small for the at-the-center location of the atom [2]. 

Another prediction made in [2] concerns the modification of the atomic Auger width 
due to the interference between a direct wave of the electron emitted in the Auger 
process and its waves scattered from the cage. This effect is identic to the one in the 
photoionization process of an endohedral atom, where the interference of the direct and 
scattered photoelectron waves can lead to the so-called 'confinement resonances' [8] in 
the photoionization cross section. In recent years a number of theoretical predictions 
have been made on the properties of the confinement resonances in various spherical 
endohedral systems (see, e.g., the review [9]). However, so far these predictions have 
not been supported experimentally. The explanation of the discrepancy between the 
predictions and experimental results was presented recently [10]. It was shown that 
the structure of confinement resonances in the photoionization cross section of an 
endohedral atom is very sensitive to the mean displacement (a) of the atom from the 
cage center. The resonances are strongly suppressed if 2(a) exceeds the photoelectron 
half-wavelength. Decrease in the amplitude of the confinement resonances with the 
increase of the doped atom displacement from the cage center was also noted in [11]. 
Taking into account similar nature of the confinement resonances in the photoionization 
and Auger decay processes we can state that the interference effects will not affect the 
Auger widths if the emitted electron wavelength satisfies the criterion formulated above. 

The experimental data on the photoionization of A@Cat are sparse due to the 
difficulty to produce sufficient amounts of purified endohedrals for the gas phase 
experiments [1, 12]. The cross sections in the region of giant atomic resonances were 
measured for Dy@C 82 [13], Ce@C 82 [14], Pr@C 82 [15] and Ce@C£ 2 [1,12,16]. 
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Apart from these the measured cross sections for photoionization of the 
endofullerene Sc3N@C^ in the photon energy range 30. . .45 eV were reported in [1]. 
These data are of the prime interest in connection with the topic of the present paper. 
Indeed, as it was shown theoretically [17,18] and also confirmed experimentally [19], 
the cross section of photoionization of an isolated ion Sc 2+ in the indicated energy 
range is dominated by a set of narrow resonances due to the excitations of the 3p 
electrons to intermediate autoionizing states. Hence, taking into account that each 
scandium atom in ScsNOCso carries a positive charge q ~ +2e [20], it is natural to 
expect the manifestation of the autoionizing resonances in the photoionization of the 
endofullerene. However, such distinct structure has not been seen in the experiment [1]. 
Instead, the presence of the endohedral molecule SC3N resulted in a single wide peak 
(of 6 eV full width at half maximum) in the photoionization curve of Sc3N@Cso- In [1] 
it was mentioned that such a dramatic modification can be due to the presence of the 
fullerene cage which could cause significant broadening of any atomic resonance features 
in photoionization. 

The intriguing experimental results of [1] has stimulated our present study on the 
dependence of the atomic non-radiative decay processes on the doped atom displacement 
from the cage center. For doing this we developed a formalism which allows one to 
determine the corrections to the atomic Auger width due to the dynamics of the fullerene 
electrons. Within the framework of the formalism, which is presented in section 2, one 
accounts for the off-the-center position of the doped atom and expresses the corrections 
in terms of the fullerene multipole dynamic polarizabilities of the second kind. On 
the basis of the developed formalism in section 3 we carry out quantitative analysis 
of the dependence of the autoionizing width on the displacement a. This is done 
for the Sc 2+ ion encaged in the fullerene Cgo- The obtained data are used further 
to calculate the photoionization spectrum for the endohedral complex Sc3N@Cg" in 
the photon energy range 30 ... 45 eV. In Appendix A we discuss different types of 
polarizabilities of the hollow objects and present the model description of the dynamic 
multipole polarizabilities of a fullerene. 

The atomic system of units is used through the paper. 

2. The formalism 

2.1. Auger width in an isolated atom 

For the sake of clarity and consistency of the notations and terminology adopted below 
in the paper, let us outline basic formalism related to the calculation of the Auger width 
in an isolated atom. 

The Auger decay of the vacancy in an isolated atom can be illustrated by the 
diagram presented in figure 1(a). The inner vacancy / is filled by the electron from 
the outer state i. Due to the interelectron Coulomb interaction v = l/\r — r'\ the 
released energy E{ — Ef is transferred to another electron from state j which becomes 
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ionized (p stands for the asymptotic momentum of the outgoing electron). The energy 
conservation implies £j + e,j — Ef + e, where e = p 2 /2. 




Figure 1. (a) Diagrammatical representation of the Auger decay in an isolated atom. 
Both transitions i — > f and j — > p occur with atomic electrons. Notation v stands for 
the (effective) Coulomb interaction. 

(b) Additional channel in the endohedral atom: the decay occurs via the virtual 
polarization of the fullerene (the particle- hole excitation — > n). V stands for the 
atom-fullerene Coulomb interaction, the filled circle denotes the effective vertex of the 
particle-hole excitation in the fullerene. 



The amplitude Af^^ p of the process from figure 1(a) is proportional to the two- 
electron Coulomb matrix element (fp\v \ji). Presenting the wavefunctions of single- 
electron states i,j, f as products of angular and radial parts, and using the standard 
multipole expansions of v and of the outgoing electron wavefunction, one derives the 
analytical expression for (the subscript 'A' stresses that the quantity refers to 

the isolated atom). In what follows, for the sake of concreteness, we focus on the dipole 
transitions % — > f and j — >■ p. As it is discussed further in the paper, the presence of 
the fullerene cage, which radius noticeably exceeds that of the endohedral atom, most 
strongly affects the dipole Auger transitions. In the dipole approximation the amplitude 
can be written as follows: 



a dip, A 




h 1 / 

-mj fi m 



Ylmip) , (1) 



where I a \ ° 1 stands for the 3j-symbol. The non-indexed quantities on the right- 



a (3 7 

hand side refer to the outgoing electron, whereas the indexed ones (the subscripts 
a = i,j,k) - to the states i,j,k. The factors S SfSi 5 SSj take into account that within 
the non-relativistic framework the spin projections s are not changed in the transitions 
i — y f and j — > p. The notation v a = (n a , l a ) stands for the set of principal and orbital 
quantum numbers, u = (e,l), and 5 U are the scattering phaseshifts. The quantity 
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Vvfv, ViVj denotes the dipole radial Coulomb matrix element: 

oo oo 

^w^H^f'o' o)(o o o)// drdr '^' (r) ^ (r ' ) ^ F "- (r)P ''' (r ' 



,(2) 

A=l 



where U lahlc ... = y/(2l a + l){2l b + l)(2/ c + 1) . . .. 

Equation (1) has been derived in a single-electron approximation. To account for 
the electron correlations one substitutes V^i;^- with the corresponding matrix element 
of the effective interaction (see, e.g., [21]). 

The width, which defines the probability (per unit time) of the Auger transition, 
one calculates starting from the following general formula: 



k mimrmj s i B f I 



(3) 



limfTtij 

"3 



where T VfV . UiV . stands for the partial width of the transition v^Uj — > vjv. 

Using (1) in (3), one derives the following expression for the partial width of the 
dipole Auger transition in the isolated atom: 



-pdip 



A 



47T 



2 



(4) 



2.2. Correction to the width of the Auger decay in the endohedral atom 



If the atom is encaged in a fullerene Cm, the atomic decay process can be strongly 
influenced due to the excitations (real or virtual) of the fullerene electrons [2,4]. 

To start with we mention, that the atomic Auger decay can occur via the additional 
channel, diagrammatical representation of which is given in figure 1(b) [2]. Here, the 
energy Uif = — Ef released in the transition i — > f is transferred to the electron j 
not directly, as in the process from figure 1(a), but via the virtual excitation of the 
fullerene. The diagrams 1(b) constitute the correction term, AAf®^, to the amplitude 
of the Auger decay. Hence, the total amplitude of the Auger decay in the encapsulated 
atom reads 

Physically, the virtual excitation implies that the Coulomb interaction between two 
atomic electrons, i and j, is modified due to the polarization of the fullerene shell. The 
polarization is dynamic, i.e. it depends on the transition energy Uif, so that AAf®^^ is 
proportional to the dynamic susceptibility of the fullerene. Thus, one can expect, that 
in those Wjj-regions where the modulus of susceptibility is large enough, the additional 
channel can modify noticeably the width of the transition ij — > fp. 

In [2] the additional channel of the dipole Auger decay was analyzed for an atom 
placed at the center of a spherically-symmetric fullerene. Assuming the atomic radius 
R a to be much smaller than the fullerene radius R, the authors expressed the amplitude 
AA^j^p Cn in terms of dynamic dipole polarizability of the fullerene. Then, carrying out 
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the order-of-magnitude estimate they found that AA^j® Cn / ' A^^ p ~ \R a / Rj < 1, 
i.e., the correction is much smaller that the Auger amplitude in the isolated atom in the 
whole range of u. 

Below in the paper we demonstrate, that the relative magnitude of the correction 
term is governed not only by u (this has been already noted in [2,4]) but also is strongly 
dependent on the position of the atom inside the cage. The parameter, which defines 
this dependence, is the ratio a/R, where a stands for the displacement of the atom from 
the cage center. Below in this section we demonstrate, that the correction term as a 
function of both cjj/ and a/R varies by orders of magnitude. This can result in the 
dominance of the additional channel over the direct Auger decay (see the case study 
presented in section 3.1). 

Let us outline the derivation of the correction term AA^^ ® Cn for the dipole Auger 
transition occurring in the atom (nucleus charge Z a , average radius R a ), located inside 
spherically-symmetric fullerene (average radius of the cage R, average width of the shell 
AR) and displaced by a from the cage center, see figure 2. 




R 



Figure 2. Illustrative picture of the atom located inside the spherical fullerene shell. 
The shell average radius and width are notated as R and AR, respectively. The atomic 
nucleus (the charge Z a ) is displaced by vector a from the center. Position vectors of 
the fullerene, e c , and the atomic, e a , electrons are notated as r and p. 



Two diagrams in figure 1(b) correspond to the following analytic expression: 

(f;n\V\0;i)(p;0\V\n;j) , (f;0\V\n;i){p;n\V\0;j) 



"S fSi®SSj 



\j i ><<\ v \v, i 'ixi>,v\ v vh Ji + \j r'w\/y"-i " i^j/ /g\ 

\_ w n0 — Uif U n0 + U>if 

The sum is carried out over the complete spectrum of the fullerene excited states (the 
integration over the continuous spectrum is implied), u n0 = e n — e is the energy of the 
virtual transition — > n. 

The quantity V = J2 a c l/\r c — a — p a \ stands for the Coulomb interaction between 
the atomic and fullerene electrons, see figure 2. Assuming the fullerene electrons to be 
located farther from the center than atomic electrons, i.e. r > |o + p|, and choosing the 



Vacancy decay in endohedral atoms: the role of non-central position of the atom 7 



z-axis along a one expresses V in terms of multipole series [22,23]: 



Here, the quantity Qi m , defined as 



Qlm({Pa}) = J^J J2Pa Y lm(Pa), (8) 

a 

is commonly termed as the operator of 2 z -pole electric moment of the atom (the sum is 
carried out over atomic electrons). The quantity 



c c 

also is related to the multipole expansion of the electrostatic field created by electron 
charge distributed in the fullerene shell (the sum is carried out over the fullerene 
electrons). To distinguish Ci m from Qi m we call the former as the operator of "interior 
2'-pole moment" of the fullerene (see Appendix A.l for the details). 

The argument a on the left-hand side of (6) is added to stress that the amplitude 
is dependent on the displacement from the center. 

Inserting (7) into (6), one carries out the intermediate algebra and arrives at the 
following expression for the amplitude of the dipole i — > f and j — > p transitions: 



1 l=L±l 



% / V^ n 2L~2 (2£ + l)! ~(C N )/ x / If 1 I 

M1M2 

(L-l I L\(L-l 1 L\( lj 1 l\ 

x n m n /i# *hn(p)- 10 

\ Hi —M J \ — /i 2 M / \—mj fi 2 m ) 

Here and dw] stand for the radial matrix elements of the atomic dipole moment, 
a^ n \ujif) is the fullerene's dynamic 2 L -pole polarizability of the second kind. This 
type of polarizability appears when one is interested in the modification of the electric 
field inside the fullerene due to its polarization under the action of the 2 i -pole external 
field whose source is also located in the fullerene interior. This is in contrast to the 
"conventional" polarizability a L N (uiif) which is responsible for the same effect but 
in the case where both the source and the observation point are located outside the 
fullerene (see Appendix A for more details). 

Summing up (1) and (10) one constructs the total amplitude A^j^f ® Cn of the dipole 
Auger decay. Then, using A^^® Cn in (3), one derives the following set of expressions, 
which defines the partial width of the transition v{Uj — > v$v in the encaged atom: 

pdip ) A<BCW \ _ pdip.A p ( \ 
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The factor F UfU . UiUj (a), which depends on the atom's displacement, takes into account 
the width modification due to the presence of the fullerene (for an isolated atom 
1). It can be written in the form 



The factor k 



1 -2k 



UfU\UiUj 



Re Tii(cUif] a) + k 



2 



( 



T, 2 (u>if, a) 



.(12) 



is constructed from the atomic characteristics only: 



K 



F7.A 

v f v i 

^dip,A 



1/2 



(13) 



where is the partial radiative width of the transition Vi — > Uf, a v . (u)if) is the partial 
cross section of photoionization of the subshell Uj, and T^)^, iVj is given by (4). 

The quantities Ei^cfj/; a) are dependent on the displacement a and on the 
fullerene's dynamic multipole polarizabilities of the second kind: 

-2 



£i(wi/; a) 
T, 2 (uif, a) 




a 



(On), 



if J 



(14) 



where c ~ 137 is the light velocity in the atomic units. As a function of the fullerene 
cage radius R the polarizability aj .\ouif) scales as ~ R 2l ~ A (see (A. 7) and (A. 12)). 
Therefore, the expansion parameter in the series (14) is (a/R) 2 . 

In the limit a = (atom at-the-center) only dipole terms with / = 1 contribute to 
the sums (14). Hence, one finds S 2 (tUj/; 0) = and 



Ei(wi/;0) 




7T U, 



2 5f N) i 

if 



"if) 




c a 



(On), 



if) 



i? 6 



(15) 



where a^- CN \uif) is the dynamic dipole polarizability of the fullerene (the transformation 



from a\ 



(On), 



>f 



to a 



(On) 



if), 



which leads to the approximate equality, is explained in 



Appendix A. 6). As a result, the factor (12) reduces to 



1 - K 



VfV\ViVj 




a 



(On) 



UJ 



if) 



if 



R 6 



(16) 



Expression (16) was obtained in [2]. The authors, estimating the magnitude of the 
correction term in the brackets as (R a / R) 3 1, derived F UfU]U . u .(0) w 1. Thus, it was 
tacitly concluded that the additional channel of the atomic Auger decay, figure 1(b), 
does not affect the decay rate if the atom is placed at the center. 

However, as the displacement a from the cage center increases, the larger 
contribution to the sums from (14) comes from the terms with higher I. This leads 
to a noticeable increase in magnitude of the factor F VfV . ViV . (a) not only in the vicinity of 
the fullerene giant plasmon resonance but in a much wider range of transition energies 
ujfi. Potentially, the atomic Auger decays with ujfi within the interval 10 ... 40 eV can 
be affected. 
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To illustrate this statement, in figure 3 we plot the dependences of /i(u;;a) 



— 2Re£i(u;;a) and / 2 (u;;a) 





2 




Si(cj;a) 


+ 


S 2 (w;a) 



/2 (see (12)) on the transition 



energy lo and for several a/R ratios. The calculations were performed for the fullerene 
ion CgQ . For the cage radius we used R = 4.15 A [24], the width of the shell was fixed 
at 1.5 A. The fullerene polarizabilities a\ CN \u) were calculated within the framework 
of the plasmon resonance approximation (see, e.g., [25,26]) and with account for two 
coupled surface plasmons excited in a fullerene (see, e.g., [27-29] and references therein). 
The details of the PRA formalism in describing the dynamic polarizabilities of the shell 
of a finite width are given in Appendix A. 4. The list of cage-dependent parameters, 
used in the calculations, one finds in Appendix A. 5. 




30 40 
transition energy (eV) 



30 40 
transition energy (eV) 



50 



Figure 3. Dependences fi(uj;a) = — 2ReEi(w;a) (left panel) and /2(w;a) = 

2 2 

£i(w;a) + £2(^5 a) /2 (right panel) on the transition energy w calculated for 
various a/R values as indicated (see equations (12) and (14)). The data refer to 
the fullerene CgQ~ (see Appendix A. 5) and explanations in the text. 



By itself large absolute values of the functions fi(u; a) and/or fi(u)\ a) in the energy 
interval ~ 10 . . .40 eV do not imply that the last two terms on the right-hand-side of 
(12) will result in a strong modification of the Auger width. Indeed, these correction 
terms are dependent also on the factor n UfV - UiVj , whose magnitude is determined by 
atomic characteristics taken at given energy (see equation (13)). Nevertheless, studying 
the curves in figure 3 one can speculate that, as the displacement increases, the region 
u « 30 . . . 40 eV becomes favourable from the viewpoint the Auger width modification 
due to the channel depictured in figure 1(b). Within the indicated interval not only 
strong inequality f 2(00; a) 3> 1 is valid but also f 2 (u;a) > 0, which means that the 
interference of the two Auger decay pathways, figures 1(a) and 1(b), is positive. This 
leads to the enhancement of the decay rate. 

For a centrally positioned atom the enhancement is more pronounced at u ~ 20 eV 
(see the solid curve in both graphs in figure 3. This energy corresponds to the dipole 
plasmon resonance associated with the symmetric mode of the coupled oscillations 
of the charges induced at the fullerene surfaces (see Appendix A. 4). The second 
surface plasmon peak, which is due to the asymmetric mode of the oscillations, is less 
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pronounced in the profile of the dipole polarizability of the second kind 5r N >) (see 
figure Al). As the multipolarity / increases, the resonance energies move towards each 
other and their intensities gradually equalize. For I ^> 1 the resonances merge at u ~ 30 
eV. This explains the increase of fi(uj;a) and / 2 (^;a) with a: as the displacement 
increases, the terms with / ^> 1 contribute more to the series (14). 

Therefore, a good candidate for the Auger width modification due to the mechanism 
from figure 1(b) would be an atom whose Auger spectrum has pronounced peaks in the 
range u ~ 30 ... 40 eV, and which, being encaged, is located far off-the-center. The 
molecule Sc3N@C^ is a good example of this sort. Indeed, three scandium atoms are 
symmetrically located far from the center being displaced by a ~ 2 A [30,31]. In 
this complex, each scandium atom donates 2 electrons to the fullerene cage. Hence, the 
theory developed above in this section can be applied to the Auger transitions in the ion 
Sc 2+ , for which the Auger spectrum has been investigated both experimentally [19] and 
theoretically [17,18,32]. In section 3.1 we consider the modification of the Auger widths 
in Sc 2+ due to the presence of the CgQ shell. The obtained results are applied further in 
section 3.2 to calculate the photoionization spectrum of the Sc3N@Cg molecule. The 
results of our calculations are compared with recent experiment [1]. 

2.3. Additional channel of the Auger decay in the endohedral atom 

In addition to the modification of the atomic Auger decay via the change of the electric 
field at the atom, the presence of the fullerene shell opens another channel of the Auger 
decay, which is absent in the isolated atom [2,4,5]. Diagrammatical representation 
of this process is given in figure 4. The energy Uif, released in the atomic transition 
is transferred, by means of the Coulomb interaction, to the fullerene electron 
which becomes ionized. In [4] the process was called interatomic Coulomb decay, thus 
pointing out its similarity with the interatomic decay in smaller systems (e.g., in neon 
dimer [6]). 




p 

Figure 4. Decay of the atomic vacancy i — > f via the fullerene ionization — > p. V 
stands for the Coulomb atom-fullerene interaction, the filled circle denotes the effective 
vertex of the particle-hole excitation in the fullerene. 

It was pointed out in [2,4] that this decay channel in the endohedral atom becomes 
especially important if the non-radiative transition % — y f is energetically forbidden in 
the isolated atom. In this case one compares the additional Auger decay width T^®^ N 
and the width of the radiative decay, T^ A . The estimates, carried out in [2,4] for the 
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atoms placed at-the-center of the Cqq cage, show the dramatic increase of the decay rate 
(by the factor r^®^ N /r^ A ~ 10 5 — 10 6 ) for the transitions with ojfi lying in vicinity of 
the giant plasmon resonance at 20 eV. 

In this section we demonstrate, that for off-the-center atoms the increase of the 
decay rate can become even more pronounced and in much wider range of the transition 
energies. 

To construct the amplitude of the process presented in figure 4, we use the 
approximate formula (7) for the atom-fullerene Coulomb interaction. The amplitude is 
used further in (3) to derive the width. Restricting ourselves to the dipole transitions 
z/j — > Uf, we obtain the following expression for the partial width of the transition as a 
function of the atom displacement: 

r^ A@CN (a)=^, /; a)r^. (17) 

The factor /F(ujif\a) stands for the ratio of the Auger-decay width in the endohedral 
atom to the radiative width in the isolated atom. It depends on the displacement and 
the transition energy, and can be written as follows: 

3 00 

T{uj if ;a) = V/(2/ + l)a 2 '- 2 Im5^(^ f ). (18) 

luj h i=i 

Similar to the series from (14), the expansion parameter is (a/R) 2 . 

Placing atom at the center and accounting for the approximate relation (A. 19), 
one derives J 7 ^/; 0) = {3c 3 /2uf f ) lma^ N (u lf ) « (3c 3 /2R 6 uf f )lma^ N {u if ). Using the 
latter expression in (17), one arrives at the formula derived in [2]. 



a/R=0 

\ a/R=0.1 




transition energy (eV) 

Figure 5. Dependence of .F(a;u;) from (18) on the transition energy oj calculated for 
several a/R values as indicated. The data refer to the fullerene Cg^ (see Appendix 
A. 5) and explanations in the text. 

Dependence of the ratio J-(oo; a) on the transition energy, calculated for several a/R 
values, is presented in figure 5. The data refer to Cgg . The parameters of the fullerene 
cage and of the dynamic polarizabilities a^ N (u) are given in Appendix A. 5. 
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It is worth noting change in the F(uj; a) profile with increase of the displacement 
a. For the at-the-center atom (or slightly off-the-center) the behaviour of ^(u); a) is 
determined by the imaginary parts of a^ N (u) with / m 1 (see lower left graph in figure 
Al). As a result, there is a pronounced maximum at ~ 20 eV, - the energy of the first 
surface plasmon excitation. In the vicinity of the second surface plasmon peak (f» 36 
eV) the value of ^(oj; a) is less by an order of magnitude. As the displacement increases, 
the higher-/ terms contribute to the right-hand side of (18). For / ^> 1 both plasmon 
modes merge, so that for each I the dependence Im a t N (u) has a single peak at u ps 30 
eV. Hence, in this energy range the ratio J 7 (a;; a) acquires additional enhancement. The 
curves with a/R = 0.4. . .0.6 show that large enhancement can be achieved in a wide 
interval of energies: u ~ 20 ... 40 eV. 

Therefore, we conclude that the decay channel, presented by figure 4, can be 
additionally intensified (by more that the order of magnitude) if the endohedral atom 
is located sufficiently far away from the cage center. 



2-4- The radiative decay in the endohedral atom 

The radiative decay i — > f + 7 in an isolated atom is represented by the first diagram in 
figure 6. The correction due to the virtual polarization of fullerene under the joint action 
of the Coulomb interaction and the photon field is illustrated by other two diagrams in 
the figure. For a centrally-positioned atom this process was considered in [2]. 




\ \ 

Figure 6. The radiative decay i — > f in an isolated atom (the first diagram) and the 
correction due to the virtual excitation — > n of the fullerene (two last diagrams). 
The dashed line depicts the emitted photon, other notations as in figure 1. The energy 
conservation law implies uj = ioi f . 



For a non-central position, applying the approximations described above, one 
derives the following formula which relates the partial widths of the (dipole) radiative 
decay in the endohedral atom and in the isolated one: 

r#? ° N = V{u if ) rz+ , (19) 

where 

(20) 



/3i° N (^/) 
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Here /3 1 N (ojjj) stands for the dipole shielding factor of the fullerene. This quantity 
determines the modification of the electric field inside the fullerene due to its polarization 
under the action of the external uniform electric field (see Appendix A for more details). 
Using the approximate relationship between /3 1 N (o; i /) and the dipole polarizability of 
the fullerene, P? N (oJ if ) « (u if ) / R 3 (see (A.19)), one reduces the right-hand side of 
(19) to the formula derived in [2] for the atom at-the-center. 

The quantity T>(u) is a so-called dynamical screening factor [27, 28, 33]. It is 
equal to the ratio of the intensity of the total electric field E to t at the atom to the 
intensity of the external field E applied to the endohedral system A@Cn- The total 
electric field E to t = Eq + AE contains additional term AE which is due to the shell 
polarization. In the cited papers it was demonstrated that T>(oj) defines the increase of 
the photoabsorption rate of the atom due to the presence of the fullerene shell. Radiative 
decay is the inverse process for the photoabsorption. Therefore, the same dynamical 
screening factor relates the radiative width in the encaged and isolated atoms. 

The factor T>(u), as it is defined by (20), does not depend on the distance a from the 
cage center. This is a consequence of the approximation made to derive equation (20). 
Namely, when calculating the additional electric field AE we ignored the reciprocal effect 
of the polarized atom on the shell. Thus, it was assumed that the shell is polarized only 
under the action of the uniform field E (u) of the dipole photon. The field E (u) gives 
rise to a uniform polarization of the field inside the cavity of the fullerene. As a result, 
the dynamical screening factor does is not sensitive to the position of the endohedral 
atom [27,33]. Such an approximation must be modified if the atomic dynamic dipole 
polarizability a^(cu) is large enough to cause a noticeable additional polarization of the 
shell. In this case, as it was demonstrated in [28,34], the dynamical screening factor 
acquires dependence on the displacement vector a. In the present paper, for the sake 
of simplicity, we ignore this reciprocal effect of the polarized atom on the shell. 



40 




photon energy (eV) 



Figure 7. Dynamical screening factor V(u)) = \l- fc N (w) for the fullerene C 
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The dependence of the dynamical screening factor on the photon energy is 
presented in figure 7. The calculations, performed for Cgg, were carried out within the 
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plasmon resonance approximation, see Appendix A. 4. The profile of this dependence, 
characterized by a powerful peak in the vicinity of the symmetric surface plasmon energy 
and with the extended right shoulder due to the activation of the antisymmetric plasmon 
excitation, is similar to those reported earlier for other spherically-symmetric fullerenes 
of finite thickness [27,34]. Thus, one can expect the order-of-magnitude increase in the 
decay rate magnitude for radiative transitions within « 20 ... 25 energy range, and less 
than that for u m 25 . . . 35 eV. 

To conclude this section let us note that the applicability of the plasmon resonance 
approximation for the dynamical screening factor was tested in [28, 34, 35] against 
the TDLDA calculations [36,37]. A good quantitative agreement between these two 
approaches was reported. 

3. Numerical results 

In this section we apply the approach described above to carry out model calculations 
of the photoionization cross section of the endohedral complex Sc3N@C^ in the photon 
energy range 30 ... 45 eV. Within this range of photon energies the photoionization 
spectrum of an isolated ion Sc 2+ is dominated by autoionizing resonances (of the widths 
less than 1 eV) due to 3p excitations [17-19]. However, recent experiments by Miiller et 
al [1] demonstrated that a distinct resonance structure is absent in the photoionization 
spectrum of the ion Sc3N@Cg" . 

In what follows, we start with the quantitative analysis of the dependence of the 
autoionizing widths in the encaged Sc 2+ on the displacement of the ion from the cage 
center due to the mechanisms discussed in sections 2.2 and 2.3. The obtained data are 
used further to calculate the photoionization spectrum of the endohedral complex. 

3.1. Auger widths in Sc 2+ and in Sc 2+ @Cgg 

Cross sections for the photoionization (PI) of Sc 2+ ions were measured by employing the 
merged ion-photon beams technique [19]. It was noted (see also [18]), that in addition to 
a direct photoionization pathway, the process can proceed via an intermediate resonance 
state according to the scheme 

hu + Sc 2+ — ► (Sc 2+ )* — ► Sc 3+ + e~ . (21) 

The intermediate state decays via (super-) Coster-Kronig transition resulting in a 
resonance line in the PI spectrum. 

The Sc 2+ ion beam, used in the experiment, contained ions in the ground state 
(the configuration [Ne]3s 2 3p G 3d 2 D 3 / 2 ) and in the first two excited, metastable states 
{[Ne]3s 2 3p 6 3d 2 D 5/2 and [Ne]3s 2 3pHs 2 S 1/2 ) as well. This lead to a variety of the 
intermediate states (Sc 2+ )* and, as a result, to a number of measured peaks in the 
spectrum within the photon energy range 30 ... 45 eV. The parameters of experimentally 
measured peaks is contained in tables and figures in [19]. Theoretical data on the 
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resonances together with the theory-versus-experiment discussion one finds in [18] (see 
also [32]). 

For our case study we restrict ourselves to the transitions from the ground state 
only. In table 1 the data are presented on most pronounced PI resonances within the 
interval uj = 31 . . .42 eV. For each excited state (Sc 2+ )* the following information is 
included: resonance energy uj 2 i (subscripts " 1" and "2" mark the ground and the excited 
states, respectively), peak value a of the PI cross section, Auger width r^ Sc2+ ^*, peak 

1 /2 

area S, radiative width and the parameter k = (t^^^* a /T 2 ^i 1 calculated in 
accordance with (13). 

The values of lo 2 i, c, r( Sc2+ )* and S were deduced from the data presented in Table 
II and Figs. 4,5,8-10 in [19]. When indicated, the Auger width (or the cross section 
a) were calculated assuming the Lorentz profile of the peak. This was done using the 
relation crr( Sc2+ )* = 2S/n. Thus, in our calculations we ignored the asymmetry seen in 
the experimentally measured peaks. 



Table 1. Transition energies u>2i from the excited state (Sc 2+ ) to the ground state 
Sc 2+ (Z?3/2), peak widths r' Sc + ) , peak values of the PI cross section a, peak areas S, 
radiative widths T^^i and parameters k (see equation (13)). 



(Sc 2+ )* a 


W21 

(eV) 


r (sc 2 +r 

(meV) 


S 

(eVMb) 


a 
(Mb) 


F 7 

(10" 6 eV) 


K 


3d" 2 F 5/2 


31.66 


116 


0.87 


4.8 b 


0.15 


4.71 x 10" 


-4 


3rf4s 2 Pi /2 


33.22 


45 


1.22 


17.3 b 


0.70 


3.10 x 10" 


-3 


3d4s 2 F 5/2 


34.73 


53 


5.40 


64.9 b 


1.13 


7.02 x 10" 


-3 


3d" 2 F 5/2 


37.14 


847 


47.9 


36.0 


11.4 


4.17 x 10" 


-3 


3d 22 P 1/2 


39.63 


11. 5 C 


10.1 


560 c 


8.24 


0.120 




3d" 2 P 3/2 


39.72 


11.3 d 


1.60 


90 c 


0.65 


1.36 x 10" 


-2 


3rf 22 D 5/2 


40.22 


2.3 d 


1.60 


450 c 


0.45 


5.60 x 10" 


-2 



a [Nc]3s 2 3p 5 is omitted in the excited-state designation. 

b Calculated from the Lorenz profile (see explanation in the text). 

c Deduced by digitalizing Figs. 9 and 10 in [19]. 



The radiative width of the transition 2 — > 1 was calculated as follows 

1^ = ^/1- (22) 

where g\ — 4 and g 2 are statistical weights of the ground and the excited states, and 
f\->2 is the oscillator strength of the transition 1 — ?- 2. For a given transition /i^ 2 is 
related to the PI peak area S through /i_> 2 = ( c /2tt 2 ) S. 

Once the values of Tj^i and k are established one can use them in (11), (12), (17) 
and (19) to calculate the total widths of the excited states for the ion Sc 2+ encaged in 
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a fullerene. The results of such calculations, carried out for the system Sc 2+ @Cg , are 
presented in table 2. 



Table 2. Total widths T tot 80 (a) for several a/R ratios. For each transition the 

upper line is due to the ^2^1 an d K data from table 1, the lower line was obtained 
accounting for the correction due to Sossah et al [18] (see explanation in the text). 



(Sc 2 +)* 


W21 
(eV) 


r (Sc 2 +)* 

(meV) 


a/R = 


r (Sc 
1 tot 

a/R = 0.4 


' 80 (a) m eV 
a/R = 0.5 a/R = 0.55 


a/R = 0.6 


3rf 22 F 5/2 


31.66 


116 


0.116 
0.115 


0.118 
0.120 


0.127 
0.135 


0.141 
0.158 


0.177 
0.223 


3d4s 2 P 1/2 


33.22 


45 


0.048 
0.049 


0.071 
0.090 


0.117 
0.175 


0.183 
0.302 


0.351 
0.650 


3d4s 2 F 5/2 


34.73 


53 


0.062 
0.069 


0.114 
0.165 


0.215 
0.370 


0.361 
0.685 


0.761 
1.595 


3rf 22 F 5/2 


37.14 


847 


1.027 
1.150 


1.343 
1.702 


1.829 
2.599 


2.432 
3.768 


3.849 
6.678 


3rf 22 Pi/ 2 


39.63 


11.5 


0.118 
0.246 


0.363 
0.837 


0.931 
2.267 


1.875 
4.685 


4.769 
12.19 


3rf 22 P 3/2 


39.72 


11.3 


0.018 
0.022 


0.028 
0.041 


0.044 
0.073 


0.066 
0.119 


0.122 
0.243 


3rf 22 D 5/2 


40.22 


2.3 


0.008 
0.013 


0.020 
0.038 


0.044 
0.094 


0.082 
0.186 


0.195 
0.468 


The first three 


columns 


in table 


2 are identical to those 


from table 


1 and are 



reproduced for the sake of convenience only. The presented values of the total width 
Ptot ^ @ ° so ( a ) were obtained by summing three terms, calculated from (11), (17) and 
(19) for several values of the a/R ratio as indicated, f 

For each excited state (Sc 2+ )* the upper line presents the values of ^ @Cso (a) 
obtained using the data on r 2 _>! and k from table 1. The lower line was obtained using 
the r 2 _,.i and k values corrected due to the recommendation made by Sossah et al [18]. 

In the cited paper theoretical calculations of the PI cross section of Sc 2+ were 
performed for photon energies from threshold to 68.0 eV. The result of calculations 
were compared to the experimental data from [19]. One of the results of this comparison 
concerns the discrepancy in the calculated and measured values of the total oscillator 
strength of the 3p 6 shell, which must be equal to 6 in accordance with the Thomas- 
Reiche-Kuhn sum rule. On the basis of physical arguments, presented in pp. 7 and 8 
of [18], the authors indicate that "... the total oscillator strength in the photoionization 
cross section of Sc 2+ from threshold to 68 eV should be a bit under 6...". They 

I The contribution of the modified radiative width, equation ((19), to T tot 80 (a) is negligent 

compared with the other two terms. It is mentioned for the sake of consistency only. 
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mentioned further, that the calculated total oscillator strength is 5.29, whereas the 
one which follows from the experimental data equals to 3.24, which is too low. The 
conclusion, which is drawn by the authors, is as follows: "Thus we believe that the 
overall magnitude of the measured cross section is too small and should be multiplied 
by a factor of 5.29/3.24=1.63 to bring the oscillator strength to a reasonable value." 

One easily verifies that such multiplication increases and k by the same factor 
of 1.63. The lower lines in table 2 were obtained using these enhanced values. 



The data in table 2 illustrate the sensitivity of the width ^ @C80 (a) both on 
the atom displacement and on the excited state. While for all excited states (Sc 2+ )* 
the widths increase with a, the increase rate is quite different for different states. It can 
be quite moderate as a/R varies from to 0.6 (e.g., the increase by approximately 50 
% in the case of 3<i 22 F 5 / 2 at uj^i = 31.66 eV) or very sharp (up to 10 3 increase for the 
3d 22 Pi /2 state). 

Let us briefly describe the relative contribution of the two mechanisms, presented 
by figures 1(b) and 4, to the total width r tot 80 (a). With the exception for the 
3<i4s 2 Pi/2 excited state, the contribution due to the diagram 1(b) dominates by the 
factor from 2 to 20 depending on the state and on the a/R value. In the case of 
3d4s 2 Pi/2 both channels contribute to the total width almost equally. 




Figure 8. Total width ^ @ c so ^ f |- ne exc ited 3d 2 2 F state (excitation energy 

W21 = 37.14 eV) versus a/R. The dashed curve was calculated using the Tj^i and k 
data from table 1. The chained curve accounts for the 1.63-correction due to Sossah 
et al [18] (see explanation in the text). Horizontal solid line marks the width in the 
isolated ion Sc 2+ . The vertical line marks the ratio a/R — 0.55 which can be estimated 
as the largest available for Sc 2+ in Sc3N@C§o (see section 3.2). 



In figure 8 the dependence T 



(Sc 5 
tot 



■)*@C|o 



is plotted for the excited 3d F state 



(excitation energy u>2i = 37.14 eV) which is responsible for the widest peak in the PI 
cross section from the ground state of Sc 2+ [19]. In the encaged ion the width can attain 
noticeably larger values than in the isolated ion (the horizontal line in the graph), and 
this influences the profile of the PI spectrum. 
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To analyze the influence of the fullerene cage on the PI spectrum of Sc 2+ one can 
use the following simple model. The cross section a Sc 2+(u; a) is calculated as the sum of 
resonance terms corresponding to the intermediate excited states (Sc 2+ ) * listed in tables 
1 and 2. Each term is approximated by a symmetric Lorentzian line profile. Thus, the 
model does not account for the asymmetry of the peaks (see the discussion in [19]). 
Hence, a Sc 2+(co; a) was written as follows: 

ht ,. (23) 



a Sc 2+(u;;aj 



2ir ^ (uj 

j=i 



r?/4 



where j enumerates the excited states. For each (Sc 2+ )* the values of resonance 
frequency cuj, peak area Sj and peak width (for the isolated ion as well as for the 
encaged one) one finds in the tables. 
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Figure 9. The PI cross section of the isolated Sc 2+ ion (thin solid line) and of the 
Sc 2+ ione encaged in Cg^~ at different distances from the center (as indicated). 



Figure 9 presents the dependence of a Sc 2+(u;a) on u calculated for the isolated 
scandium ion (thin solid curve) and for Sc 2+ encaged in Cgg at different distances 
from the center. To calculate the resonance terms in (23) we used, for each indicated 
value of the ratio a/R, T^l ^ @Cso (a) listed in the upper lines in table 2. Within the 
model framework we assumed the peak areas S (see the third column in table 1) to be 
independent on a. Also, we disregarded possible redistribution of the ionic oscillator 
strength due to the presence of the fullerene shell [12]. 

The focus of our study is on the modification of the PI spectrum profile due to 
the Auger widths broadening. From this viewpoint the comparison of different curves 
in figure 9 is quite illustrative. For moderate values of a/R, similar to the case of 
the isolated ion, the spectrum is dominated by well-separated lines. As a increases, 
the lines, loosing in the peak values and gaining in the widths, become less accented. 
For a/R = 0.5 . . . 0.6 (this range of the a/R ratio is quite realistic, see section 3.2) the 
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spectrum noticeably flattens. This effect is most pronounced for u < 35 ... 37 eV, where 
the effect of the fullerene shell polarization on the ionic Auger decay is very strong (see 
figures 3 and 5). 

3.2. Photoionization spectrum o/Sc3N@C^q 

In [1] the experimental results for the PI cross section were reported for Sc3N@C^ . 
To carry out theoretical investigation of the PI process of this endohedral system we 
assume that both its electronic structure and isomeric form are close to those of the 
neutral compound Sc3N@Cso- which is a member of a family of the tri-metallic-nitride 
endohedral system [38]. These objects are rather curious in that often both of the two 
subsystems are unstable in isolation. When these two subsystems are combined into an 
endohedral fullerene, there is a series of electron transfers between various components 
of the subsystem that results in a mutual stabilizing effect [39] . 

Electronic structure of Sc3N@Cso metallofullerene was discussed in [20]. Each 
scandium atom donates 2 electrons to the fullerene cage. Additionally, a partial 
charge of 0.4 is donated to the nitrogen atom. The endohedral complex is therefore 
(Sc 2 ' 4+ ) 3 N 1 - 2+ @CgQ . For the purposes of the present paper we ignore the partial charges 
in all intermediate stages of the calculations. Hence, each scandium ion is treated as 
doubly ionized, Sc 2+ , and the nitrogen is considered as a neutral atom. At the final 
stage, when constructing the PI cross section of the compound Sc 3 N@Cg the partial 
charges will be accounted for following the phenomenological arguments presented in [1] . 

The isomeric form of the CgQ cage becomes of the lh type, i.e. 'nearly spherical'. 
Basing on the study of Nakao et al [24] (see also [40]) one can use the value R = 4.15 
A for the average radius of the isomer. 

All four atoms of the SC3N complex lay in a single plane [30,31]. The nitrogen atom 
forms the central part of the structure, and the three scandium atoms are positioned at 
an average separation of 2.0 A. Therefore, had the nitrogen atom resided at the cage 
center the (average) displacement of each Sc 2+ would be a = 2 A, i.e., a/R ~ 0.5. 
However, analyzing the NMR spectra of Sc3N@Cgo, Stevenson et al [31] indicated that 
the complex is not localized in any particular site in the fullerene. Therefore, it is 
meaningful to estimate the largest displacement a max which can be experienced by Sc 2+ 
in the C§q cage. Calculating the density of electron cloud in the ground state of Sc 2+ 
(for example, by means of the Hartree-Fock code [41]), one finds that -Rsc 2 + ~ 1-1 A is 
a good estimate for the radius of the ion. Taking into account that the thickness of the 
fullerene electron cloud is 1.5 A [42] for Ceo, we find the shortest distance between Sc 2+ 
and the cage as 1.1 + 0.75 = 1.85 A. This value suggests that the largest displacement 
of Sc 2+ from the center is a max « R — 1.85 A = 2.3 A, yielding a max /i? ~ 0.55. 

Hence, to estimate the increase in the Auger widths in Sc 2+ due to the excitation 
of the plasmon oscillations in the fullerene shell one can use the a/R ration from the 
range 0.5 . . . 0.55 (see table 2 and figures 8 and 9). 

In the photon energy range 30 ... 45 eV the PI cross sections of the fullerene cage [1] 
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and of the nitrogen atom [43] are smooth functions of u. Therefore, one would expect 
that any peculiarity in the cross section o" Sc3N @ c + of the endohedral complex is due to 
the scandium ions. Ignoring the difference in the PI cross section between the C^ and 
CgQ cages§ one can write the following approximate relation: 

Aa(uj; a) = <r Sc3NOC + - <r c + (w) « cr N (u) + 3a Sc 2+ (w; a) . (24) 

The experimental data on c Sc . iN @ c + and <r c + (u) one finds in figure 3 in [1]. Therefore, 
to carry out theory-versus-experiment comparison it is necessary to evaluate the right- 
hand side of (24). To calculate the term er N (cj) we used the tabulated data [44]. The 
calculation of a Sc 2+(u;a) was carried out in accordance with (23). Additionally, the 
obtained a Sc 2+(u;a) curves were shifted on the photon energy axis by +1.5 eV. This 
shift accounts (approximately) for the increase in 3p subshell binding energy of the Sc 
charge state 2.4+ instead of 2+ [1]. 




Photon energy (eV) 



Figure 10. Theoretical versus experimental data on Aa(uj;a) (see equation (24)). 
Closed circles with error bars: digitalized data from figure 3 in [1]. Thin solid curve: 
cross section Act(w; a) with the term 3<7g c 2+(cj; a) calculated for the isolated Sc 2+ ion. 
Thick curves: calculations carried out for Sc + ions displaced by a = 0.55R = 2.3 A 
from the cage center; the solid curve obtained without the 1.63-correction due to [18], 
the dashed - the correction was accounted for. See explanation in the text. 

In figure 10 we present the calculated dependences Aa(u; a) together with the 
experimental data which were obtained by digitalizing the upper graph from figure 3 
in [1]. The fitting, carried out in the cited paper, suggests that the experimental data 
can be described by a single broad peak of 6 eV full width at half maximum (this profile 
is not drawn in figure 10). This behaviour clearly differs from that of the thin solid 

§ Such an approximation seems reasonable, since 6 extra valence electrons do not affect noticeably the 
dynamics of the surface plasmons which define the magnitude and behaviour of the cross section in the 
indicated w-interval. 
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curve, which was obtained by using the PI cross sections for the isolated Sc 2+ ion in 
(24). 

The correspondence between theory and experiment can be considerably improved 
by accounting for the shell polarization and placing the scandium ions away from 
the cage center. Two thick curves in figure 10 represent Aa(u; a) calculated for 
a = 0.55i? = 2.3 A. The solid curve was obtained by using the S values from table 
1 and the upper values of ^ @C *° (a) from table 2. The dashed curve was obtained 
by accounting for the correction due to Sossah et al [18]. Within this scheme the peak 
areas S were multiplied by 1.63 and the lower values from the a/R = 0.55 column in 
table 2 were used. 

Taking into account the approximate framework used to describe the dynamic 
polarizational effects in the fullerene shell (the resonance plasmon approximation 
accompanied by several parameters introduced "by hand") as well as several simplifying 
and phenomenological assumptions made to calculate the right-hand side of (24), we 
would not claim that these curves reproduce in full the experimental data. The main 
drawback of the calculated profiles is that they still contain several peaks associated 
with the ionic transitions indicated in table 1. These features has not been seen in the 
experiment [45] . Nevertheless, the general trend, which was intuitively formulated in [1] , 
is quite clear: the presence of the fullerene cage could result in significant broadening 
of atomic/ionic resonances in photoionization. The case study, which has been carried 
in this paper and which is illustrated by figure 10, indicates that the impact of the 
cage polarization on the PI peaks strongly depends on (a) the location of the target 
atom inside the cage, (b) the energy of the atomic transition, (c) the magnitudes of the 
radiative and non-radiative width of the transition. 

4. Summary 

In summary, we have demonstrated that the Auger decay rate in an endohedral atom 
is very sensitive to the atom's location in the fullerene cage as well as to the energy u> 
released in the decay process. 

Two additional decay channels have been considered, which appear in an endohedral 
system and lead to the modification of the rate. Firstly, there is a correction to the 
direct Auger-decay channel due to the change in the electric field at the atom caused by 
dynamic polarization of the fullerene electron shell by the Coulomb field of the decaying 
vacancy. In principle, this correction can be of either sign (depending on u), and, thus, 
can either increase or decrease the decay rate. Within the second channel the released 
energy is transferred to the fullerene electron via the Coulomb interaction. 

The contributions of both additional mechanisms are expressed in terms of 
the fullerene dynamic multipole polarizabilities of the second kind. This type of 
polarizabilities appears when one is interested in the modification of the electric field 
in the interior of any hollow object (e.g., a fullerene) due to its polarization under the 
action of the external field whose source is also located inside the object. This is in 
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contrast to the "conventional" polarizability which is responsible for the same effect but 
in the case, when both the source and the observation point are located in the exterior. 
Thus, one can expect the modification of the decay rate to be most pronounced in those 
w-regions where the polarizabilities of the second kind are large enough. 

The relative magnitudes of the correction terms are governed not only by the 
transition energy but also are strongly dependent on the position of the doped atom. 
(This feature is absent in the radiative decay process.) The parameter, which defines this 
dependence, is the ratio of the off-the-center atomic displacement a from the fullerene 
center to the average radius R of the cage. As a/R increases the enhancement of the 
width occurs for the transitions whose energies are in the vicinity of the fullerene surface 
plasmons energies of high multipolarity. 

We have demonstrated that the correction terms as functions of both u and a/R 
vary by orders of magnitude. This can result in the dominance of the additional channels 
over the direct Auger decay and can lead to pronounced broadening of the atomic 
emission lines. 

On the basis of the developed formalism we carried out a case study for the system 
Sc 2+ @CgQ • It has been demonstrated that narrow resonances in the spectrum of an 
isolated Sc 2+ in the energy range u = 30 ... 45 eV are noticeably broadened if the ion 
is located strongly off-the-center (a/R ~ 0.5). Our model allowed us to carry out the 
quantitative analysis of the photoionization of the endohedral complex ScsNQCgo- We 
have demonstrated that due to the the non-central position of the three scandium atoms 
as well as to the fact that the multipolar surface plasmons energies of the fullerene lie 
within the indicated range, the additional decay channels can be responsible, at least 
partly, for the strong modification of the photoionization spectrum profile detected 
experimentally. The analysis of another reason of the lines broadening, which is also 
due to the presence of the fullerene shell but different from those discussed in this paper, 
will be published elsewhere. 
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Appendix A. Multipole polarizabilities of hollow objects 

In this Appendix we present general definitions and discuss basic properties of different 
types of multipole polarizabilities which characterize the response of a 'hollow object' 
to an external electric field. By a hollow object we will understand any object which 
has at least one cavity in its interior. A fullerene Cn, whose (average) radius R greatly 
exceeds the radius of valence electrons in a carbon atom, can be treated as a hollow 
object if one neglects the probability to observe electrons in most part of the fullerene 
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interior. 

Appendix A. 1. Definitions 

Any system of electric charges becomes polarized if exposed to an external electric 
field. The polarization of this object leads to an additional electric field whose potential 
(or/and strength) depend on the multipolarity of the external field, on the distance 
from the system to the observation point, and on the (multipole) polarizability, i.e. the 
quantity which describes the system's response to the external field. 

If an object has no internal cavity (e.g., an atom, a simple molecule), it is natural to 
assume that a source of the external field (e.g., a point-like charge) is located outside the 
object. Then one can distinguishes two different cases: the observation point is located 
(i) outside or (ii) inside the object. The different quantities are responsible for the 
system's response in these two cases. Conventionally (see, e.g., [46]), the characteristic 
associated with case (i) is called a polarizability, whilst case (ii) is described by the 
quantity termed " a shielding factor" . 

As discussed below, the response of a hollow object to the action of an external 
electric field can be described in terms of three different quantities. In addition to the 
two mentioned ones, the third quantity corresponds to the situation when both the 
external charge and the observation point are located inside the object. 

Let us briefly discuss the quantum formalism which one can apply for a quantitative 
description of these quantities. For the sake of simplicity, we consider the following 
model hollow object (called below 'a fullerene'): a system of N e electrons located in the 
vicinity of the spherically-symmetric ionic cage of the radius R. Let the fullerene be 
exposed to the field of a point charge q located in some point with the position vector 

Expanding the operator of the Coulomb interaction V = — J2j=i l/\ r j ~ r q\ 
between q and the fullerene electrons in spherical harmonics, one presents it in the form 
V = — q J2im 'filming, {rj}). The multipole potentials <f>i m (r q , {rj}) depends on whether 
the charge q is located in the fullerene exterior (r q ^> (rj)) or interior (r q <C {fj))'- 



The first one of these quantities, Qi m {{ r j})i is expressed in terms of so-called regular 
solid harmonics r l Yi m (r) (see, e.g. [47,48]), and is commonly associated with the 2 l - 
pole moment of a system of point unit charges (in quantum-mechanical terms, it is the 
operator of the 2 / -pole moment). This quantity allows one to calculate the multipole 




where 




(A.2) 
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potential created by the system in its exterior. The second quantity, Ci m ({rj}), written 
in terms of irregular solid harmonics r~ l ~ l Yi m (r) [48], determines the multipole potential 
in the system interior (see [47]). To avoid ambiguity we call this quantity the 'interior 
2'-pole moment'. 

Under the action of the multipolar field — q4>im{fq, { r j}) the electron cloud in the 
fullerene becomes polarized. This polarization causes the change in the potential, 
S(j>im(f), created by the electrons. Applying standard perturbation theory, one finds 
that 5(pi m (r) is given by the following second-order matrix element: 



*E 
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'Ira 



><?> {rj}) 
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lm 



(r, fa}) 
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(A.3) 



Here, (f>i m (r, {rj}) stands for the (operator of) multipole potential created by the 
electrons in the point r. Its explicit form follows from (A.l) by means of the substitution 
r q — y r. The sum is carried out over the complete spectrum of the fullerene's excited 
states n (the ground state is excluded) including the excitations into continuum, 
w„o = £n — £o is the excitation energy. 

Explicit form of the expression on the right-hand side of (A.3) depends not only 
on the location of the external charge q with respect to (rj), but also on the location of 
the observation point r. Considering different cases, one derives: 
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(A.4) 



The quantities ai, fa and aii are given by 
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(A.5) 



The first two of these are well-known in atomic physics (see, e.g., [46]). The (static) 
polarizability ai defines the change in the field, created by a system, in the case when 
both the source charge and the observation point are located in the system's exterior. 
The quantity 0i is called a shielding factor, and it defines the change in the field in 
the system's interior due to the charge re-distribution under the action of q located in 
the exterior. Atomic shielding factors appear when one analyzes the change of the field 
at the nucleus due to the polarization of the outer shells electrons under the action of 
external electric field. 

It is the third term, 5;, which additionally appears when a hollow system is exposed 
to the external field. In this case, both the source charge q and the observation point 
can be located in the cavity. To distinguish between a>i and ai we call the latter quantity 
the polarizability of the second kind. 
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The polarizability of the second kind naturally appears in the processes which 
involve the interaction of a subsystem /embedded by a hollow subsystem II. An example 
would be an endohedral complex A@Cn, where A stands for an atom. It can be shown 
that the dipole polarizability of the second kind 5i determines the van der Waals 
interaction between the fullerene and the centrally placed atom [49,50]. More general 
treatment of the van der Waals interaction between CV and the atom, arbitrary placed 
in the fullerene interior, involves the polarizabilities 5; of higher multipolarities [51]. 
Dynamical screening of an endohedral atom (see [27,33]) is another example of the 
process in which the polarizabilities 5; (more exactly, the dynamic polarizabilities 
see (A. 6)) manifest themselves. In this process the presence of the fullerene dynamically 
screens the confined atom from an external electromagnetic field, so that the atom 
experiences a field that is either enhanced or suppressed depending on the field frequency 
oj. It can be shown (see [28]) that to calculate the change in the atomic photoionization 
cross section due to the dynamical screening one has to calculate the shielding factor of 
the fullerene and the polarizabilities ai(u). Equation (14) in the main text demonstrates, 
that fullerene dynamic polarizabilities of the second kind determine the correction to 
the Auger decay rate of the endohedral atom. 



Appendix A. 2. Dynamic multipole polarizabilities 



The formalism outlined above can be adjusted to the case of the time-dependent external 
multipolar field. To do this one carries out the Fourier transform of the external field and 
then treats each Fourier component separately. As a result, each of the quantities ai, 
on and fii acquires the dependence on u (the parameter of the Fourier transform or the 
'frequency' of the external field). Explicit expressions for the dynamic polarizabilities 
of all types follows: 
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In the static limit u = these formulae reduce to (A. 5). 

Let us evaluate the asymptotic behaviour of oii(u), ai(u) and (3i(cu) in the region 
u 3> Ik {Ik stands for the ionization potential of the most bound electrons in the 
system). As a first step in deriving the asymptotic expressions one neglects cu^ in the 
denominators. The resulting sums over the complete spectrum of the system are carried 
out using standard methods of evaluation of the sum rules (see, e.g., [52]). To calculate 
the gradients of regular r l Yi m (r) and irregular r~ l ~ l Yi m (r) solid harmonic one can use 
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the formulae from Ch. 7.3.6 in [22]. The final result reads: 
lim a®(u) = -— 2 V2 Wb o \{0\Qi„<- 
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lim pfa) = --= y22cu n0 (0\Q lm \n)(n\C? m \0) = 



where r k denotes the mean value of r k in the ground state. 

The first formulae from (A.7) coincides with the one obtained earlier in [25]. In the 
dipole case, I — 1, it reproduces the well-known result, related to the Thomas- Reiche- 
Kuhn dipole sum rule. 

Appendix A. 3. Multipole polarizabilities of a dielectric spherical shell 

For reference, we present explicit expressions for the polarizabilities ati, on and /3j of a 
model hollow object: a spherical shell made of a material with dielectric function e. The 
inner and outer radii of the shell are notated as R± and R2, respectively. 

To find the polarizability a\ and the shielding factor 0i one places the external 
charge q in the shell exterior (r q > R2) and determines the additional potential S(p(r) 
due to the polarization of the shell. Multipolar components 54>i m (r) of this potential 
are proportional either to ai (if r > R 2 ) or to f3i (if r < Ri). To find ai one places the 
charge in the shell interior r q < Ri and looks for 8<pi m (r) in the inner region r < Ri . 
As a result, the one derives 
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The expression for the multipole polarizability presented in (A.8), coincides with the 
formula obtained earlier in [53]. 



Appendix A. 4- Multipole polarizabilities of a fullerene within plasmon resonance 
approximation 

Photoionization experiments on Ceo an d its positive ions reveal two giant resonances in 
their photoionization spectra in the photon energy region 20 ... 40 eV [54,55]. These 
resonances are of a collective nature and are due to the excitations of plasmons by the 
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photon field [29,55,56]. The existence of two resonances can be explained [29] in terms of 
two coupled surface plasmons excited in a fullerene when it is considered as a spherical 
layer of a finite thickness [27, 28, 34, 53, 57, 58] . 

In this section we apply this model to express each of the quantities ai and 5/ 
from (A. 8) as a sum of two resonance terms corresponding to the two surface plasmons 
modes. We call this approach as plasmon resonance approximation in accordance with 
the term used in the theory of plasmon resonances in metallic clusters (see, e.g., [25,26] 
and references therein). 

Let the fullerene valence electrons be distributed within a spherical shell whose 
inner and outer radii are R\ = R — AR/2 and R2 = R + A/2, with R standing for the 
(average) radius of the ionic core and A - for the width of the electron cloud (see figure 
2). Within the framework of PRA one can model the dielectric function e = e(u) of the 
fullerene electrons by a Lorentz-type dielectric function (see, e.g., [58]): 

U2p (A. 10) 



u 2 — Uq + ircu ' 

where T is a damping constant, and the term Uq is included to go beyond the free- 
electron gas approximation implied by the Drude model. The plasma frequency u p is 
found from Loi = 4irN e /V, where V = 47r(i?2 — ^i)/3 is the volume of the shell. 

Using (A. 10) in (A. 8), one represents the dynamic polarizabilities a>i{uj) and ai(cu) 
as follows: 

V ' „ . (A.ll) 

= (l + l)Mu)R- 21 -*, 




The factors R 2l ~ 2 and i?~ 2 '~ 4 , which stand for the mean (21 — 2)th and (—21 — 4)th 
powers of the radius of the electron cloud (cf., equation (A. 7)), conveniently expose the 
dependence of ai and 5/ on the shell radii: 

3 i-e +i M 3 1 i-e i+i 



R2l ~ 2 ~2i + i i-e ' R ~ 2l ~ 4 ~ 21 + 1 rir 2 ^ i-e ' (A,12) 

The (^-dependence of the polarizabilities is concentrated in the factors Ai(co) and 
Ai(u), which have similar resonance structure: 

Nu N 2l 



Ai(u) 
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~ ~ A.13 

N 1{ N 2l 



+ u 2 - iTuuj uj\ x + Uq - u 2 - 1T21OJ 
Here ouu and OJ21 are the surface-plasmon frequencies for a given multipolarity /. These 
pairs of eigenmodes arise from the fullerene having a finite thickness and therefore two 
surface charge densities. The eigenfrequency uu characterizes the symmetric mode of 
coupled oscillations, in which the two surface charge densities oscillate in phase. The 
eigenfrequency u 2 i (which is larger than uu) corresponds to the antisymmetric mode, 
in which the oscillations are in antiphase [28,53]. These frequencies are given by 



Vacancy decay in endohedral atoms: the role of non-central position of the atom 28 



where pi = yjl + 4/(7 + 1)£ 2 ' +1 . Note that for each multipolarity the frequencies uu 
and U21 satisfy the relation bj\ { + oo^i = °°p: which is a consequence of a more general 
sum rule for the magnitudes of surface plasmon frequencies [59]. 

The quantities Nji and Nji (j = 1, 2) are the (multipolar) oscillator strengths 
associated with the modes uu and U21. They satisfy the rule N11 + N21 = N11 + N21 = N e , 
and, thus, represent the numbers of valence electrons participating in the symmetric and 
antisymmetric oscillatory modes. These numbers are given by 

N ll = N2i = N e P -^, N2l = n u = N e ^^ . (A.15) 
2pi 2pi 

Comparing (A. 13) and (A.15) one notices, that the functions Ai(uS) and Ai(oj) differ 
only in the number of oscillators in the two modes: the expression for Ai(u) can be 
obtained from that for Ai{oS) by means of the substitution Nu -H- N 2 i- 

The quantities Tji (j = 1, 2) in (A. 13) denote the half- widths of the surface plasmons 
multipole excitations. They determine the decay rate from collective excitation mode to 
the incoherent sum of single-electron excitations. Within the framework of PRA they 
(as well as the quantity cj ) are treated as parameters which can be either deduced from 
the experimental data or calculated separately [26,60]. 

Within the PRA the shielding factor fiiiuj) (see the last line in (A. 8)) acquires the 
form [28,34]: 

W'- Ti'V-? ^'"'- (A - 16) 

The factor Bi(cu) has the following resonance structure: 

BiK u ) = 2 1 2 2 — ~ t — 2 2 — — • A - 17 

UJ U + UJq — U z — ll uU + OJq — 0O z — ll 21^ 

Note that, in contrast to Ai(u) and Ai(co), the "oscillator strengths" associated in Bi(uj) 
with the symmetric and asymmetric modes are the same and equal to Nu — N 2 i- 

Appendix A. 5. Dynamic multipole polarizabilities ai(u) and 5/ (a;) of C^q 

In this section we apply the plasmon resonance approximation to calculate dynamic 
polarizabilities cti(cu) and cti(u)) for the fullerene ion CgQ . 

Taking into account the study of Nakao et al [24] we assume spherical symmetry 
for the ionic cage and use the value R = 4.15 A for the cage radius. The thickness AR 
of the fullerene is set to 1.5 A, which was obtained in [42] for Cqq. These values of R 
and AR, accompanied by the number of valence electrons iV e = 326, lead to u p = 37.0 
eV for the plasma frequency and to £ = 0.69 for the inner to outer radii ratio. 

The value of u was chosen so as to equate the frequency of the symmetric 
dipole resonance (i.e., the resonance frequency of the first term in A\(uj), see (A. 13)) 
to the experimentally measured position u cxp of the plasmon resonance peak in the 
photoionization spectrum. Thus, the value uq = 13.9 eV was calculated from 
00^ + = w cxp , where Uu one calculates from (A. 14), and we took oo cxp = 21 eV, 
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which corresponds to the experimental maximum of the photoionization cross section 
for C+ [1,61]. 

To calculate the half-widths the following rule was adopted: Tji = ^jUjji (j = 1,2). 
For symmetric surface plasmon modes of all multipolarities the value of 71 was fixed 
at 0.16, which equals to the ratio of experimentally measured half-width r exp to the 
position u cxp of the dipole resonance [1,61]. For the antisymmetric modes we used 
the value 72 = 0.22, which corresponds to the r exp /w CX p ratio measured for the second 
plasmon resonance in photoionization of Cgd" ions [55]. 

The calculated dependences of real and imaginary parts of Ai(oj) and Ai(u>) on the 
photon energy u for several values of I are presented in figure Al. 
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Figure Al. Imaginary and real parts of Ai(uS) (upper graphs) and Ai(u>) (lower 
graphs) versus w calculated for CgQ and for several values of I as indicated. See 
explanation in the text for the parameters used in the calculations. 



We mention the following two features of the presented dependences. Firstly, as it 
is seen from equation (A. 14), with the growth of I the frequencies uu and U21 approach 
the same limit: uu goes to u p /\/2 from below whereas U21 - from above. As a result, 
the two resonances in both Alius) and Ai (u) merge as / increases. Secondly, the second 
(antisymmetric) surface plasmon peak is quite pronounced for Ai(u) (i.e., in the dynamic 
polar izability of the second kind) whereas in the case of Ai(u) its magnitude relative 
to the symmetric peak is much smaller. As I increases (and until the peaks merge) 
the heights of symmetric and antisymmetric resonances nearly equalize in the case of 
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Alius), whereas for Ai(u>) the antisymmetric peak quickly disappears. The reason for 
such difference in the resonance curves behaviour is related to the different number 
of valence electrons participating in symmetric and antisymmetric modes, see (A. 15). 
The ratio N 2 i/Nu rapidly decreases with / (the ratio equals to 0.3 for I = 1 and to 
0.04 for I = 10). Hence, the number of oscillators associated with the antisymmetric 
mode in Ai(u), which is N 2 u rapidly decreases, thus destroying the peak at u 2 i- On 
the contrary, in Ai(u) where the antisymmetric mode is defined by Nu oscillators, the 
relative intensity of its peak increases. 

Appendix A. 6. Approximate relationship between a>i(uj) and ai(u) for a fullerene 

To conclude Appendix let us mention a transformation which allows one to 
approximately relate the polarizabilities ai(u) and ai(u). To simplify the algebra we 
consider the case of a spherically symmetric fullerene with the cage radius R and the 
thickness AR. 

Assuming the strong inequality AR <C R one reduces the interior 2'-pole moment 
Qim({fj}) to regular solid harmonics r l Yi m (r) Qim{{fj}) (see equation (A. 2)) as follows: 

ClM r i} ) = ^ E ^ » J^J^h E tfU-,) = ^P(A.18) 

3 3 j 

Using this result in (A. 6) one arrives at 

<*M™W&> (A.19) 

Sometimes such approximation is used to estimate the dynamic response of a fullerene 
to the external field (see, e.g., ( [2])). However, the accuracy of the approximation is 
subject to (comparatively) low /-values and to a strong inequality A_R R. Thus, if one 
considers, for example, the fullerene Cqq, for which AR/R w 0.4, then equation (A.19) 
can be used only for the order-of- magnitude estimations and for I « 1. The accuracy will 
increase for giant fullerenes (e.g., for C960 where which AR/R m 0.1). Additionally, for 
low /-values this approximation neglects the difference between the profiles of symmetric 
and antisymmetric surface plasmon modes in 5/(w), /3/(w) and ai(u) (to compare the 
profiles in a; (a;) and in 5j(a>) - see figure Al). 
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